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Abstract 

We explore the dark matter detection prospects in the Minimal Supersymmetric Standard 
Model in the scenario where the scalar partners of the fermions and the Higgs particles (except 
for the Standard-Model-Hke one) are assumed to be very heavy and are removed from the low- 
energy spectrum. We analyse the neutrahno LSP (x?) in scenarios where the gaugino mass 
parameters are universal at the GUT scale and also the case where they are non-universal. This 
analysis is carried out in the framework of a Xenon-like 100 kg experiment. In general, an 
important fraction of the parameter space giving rise to the dark matter relic density measured 
by WMAP can be probed and excluded in the case of not detecting any WIMP. In the opposite 
case, once a WIMP signal has been found, we show that for a light Xi which is a higgsino- 
gaugino mixture it is possible to reconstruct efficiently the mass and the scattering cross-section 
of the neutralino LSP. Moreover, we show that it is also feasible to put strong constraints over 
some of the parameters of the Lagrangian, e.g. the higgsino and the gaugino mass parameters. 



1 Introduction 



There exists strong evidence that a large fraction of the Universe is dark and non-baryonic [HOE]. 
Weakly Interacting Massive Particles (WIMPs), with masses lying from the GeV to the TeV scale, 
are the leading cold Dark Matter (DM) candidates. The direct detection of WIMPs could not only 
directly confirm the existence of dark matter but would also probe the fundamental parameters 
of the underlying theory. Constraints on, or measurements of, the WIMP mass and scattering 
cross-section will be complementary to the information derived from different experiments such as 
colliders [HEIE], neutrino detectors [7], dark matter indirect detection |8j or even other dark matter 
direct detection experiments [H [TOl [HI [12] . It is therefore pertinent to examine the accuracy with 
which DM direct detection experiments will be able to measure the WIMP properties, if they are 
detected. 

Low-energy supersymmetry (SUSY) with i?-parity, and in particular the Minimal Supersym- 
metric Standard Model (MSSM), provides several well-motivated WIMP candidates which take the 
form of the Lightest Supersymmetric Particle (LSP) (see e.g. reference [T]). However, no SUSY 
particle has been found at particle colliders and no clear indication of SUSY has emerged from any 
of the current measurements. In this way, the initial theoretical motivations for SUSY, although 
still very appealing, are being questioned and the scale of new physics beyond the Standard Model 
(SM) is slowly, but steadily, drifting above the weak scale. In this vein, a SUSY framework named 
Split Supersymmetry |13l [T4l 115] has emerged, proposing a generic realisation of SUSY, where the 
scalar superpartners of SM fermions are extremely heavy, with a mass that in principle can reach 
the GUT scale. Moreover, the fermionic superpartners of the Higgs and gauge bosons could remain 
near the electroweak scale, protected by symmetries. Let us note that this type of models may 
appear naturally in some string inspired models with F-term uplifting [HI [171 [18] . In Split SUSY, 
two major features of the MSSM are maintained: on the one hand the unification of gauge couplings 
works essentially in the same way as in the MSSM. Actually, it has been pointed out [T3j[14] that, 
in order for SUSY to provide solutions to the unification problem, only gauginos and higgsinos, the 
fermionic superpartners of the gauge and Higgs bosons, need to be relatively light. On the other 
hand, in Split SUSY the lightest neutralino is usually the LSP, with a mass lying from few dozen 
GeV to the TeV scale, hence being a promising candidate for the cold dark matter. It is clear that 
in this limit where scalars are ultra heavy, the fine-tuning problem is reintroduced in the theory. 

However, it is not mandatory to have scalar superpartners at such a high scale. Indeed, from a 
scalar mass Ms of the order of ~ 10^ GeV the scalar particles are decoupled from the low-energy 
spectrum, and hence their interference with low-energy phenomenology is marginal. In this case, 
and as compared to the MSSM, just a slightly large amount of fine-tuning is needed. Let us recall 
that there is no compelling criterion to define the maximal acceptable amount of fine tuning |19| . 
and the choice of the upper bound on the SUSY scale Ms is somewhat subjective. An interesting 
feature of such a MSSM with heavy scalars [20], is that it is much more predictive than the general 
MSSM, since a large number of the free parameters comes from the scalar sector. At the same time, 
both the gauge couplings unification and the solution to the dark matter problem are maintained. 
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The present article is devoted, on the one hand, to the analysis of the neutralino-nuclei scatter- 
ing cross-section and to the study of prospects of neutralino dark matter direct detection within the 
MSSM with heavy scalars. Some previous studies [211 [23] have already been devoted to these 
subjects considering low-energy model-independent frameworks. Instead of following the aforemen- 
tioned approach, we will study the general model using different patterns of soft SUSY-breaking 
gaugino masses. In this way, we will consider different SUSY models in which the boundary condi- 
tions for the gaugino masses at the GUT scale are different from those of the universal scenario. 
On the other hand, we will study the possibilities offered by direct detection experiments to recon- 
struct: i) the mass and the scattering cross-section [U HOl [24] and ii) the fundamental parameters 
of the Lagrangian. The present paper aims to provide a complementary picture to that use in 
references [25l[26] and [27], in which the reconstruction of the MSSM with heavy scalars is studied 
in the framework of high-energy colliders, both at LHC and at ILC. 

The rest of the paper is organized as follows. In the next section we briefly review the MSSM 
with heavy scalars and discuss the neutralino dark matter within this model. In section [Sj we 
review the spin-independent direct detection techniques, focusing on the event rate and the WIMP- 
nucleon scattering cross-section for a Xenon- like experiment, in a microscopically model independent 
approach. This analysis lays the basis for the following study, in the framework of the MSSM with 
heavy scalars. In section [4] we present some results of dark matter direct detection prospects for a 
Xenon-like experiment. Section [5] is dedicated to the possibilities of reconstructing both physical 
observables and fundamental parameters of the model. Before closing, conclusions are given in 
section [H 



2 The MSSM with heavy scalars 
2.1 Definition of the model 

In the MSSM with heavy scalars, the scalar superpartners of the leptons and quarks as well as most 
of the Higgs bosons are assumed to be very heavy, at a common mass scale Ms ^ 1 TeV. The 
low-energy theory contains, besides the SM-like Higgs boson only the two higgsinos Hd) and 

the three gauginos: the bino (B), the wino (W) and the gluino (g). Omitting the gauge-invariant 
kinetic terms, as well as non-renormalizable operators suppressed by powers of the heavy scale Ms, 
the Lagrangian of the effective theory reads [IS, [H] 



hfjqjUieH* + hfjqjdiH + hfjijaH 



+ f~9^~9^ + fw^W- + ^BB + ,HleH, 
+ ^ [gud'^W'' + gis) Hu + ^ [-~9d'j''W'' + 9dB ) + h.c. 



fll 



where cr" are the Pauli matrices, e = i cr^ and i, j are generation indices. The SM-like Higgs doublet 
is a linear combination of the two MSSM Higgs doublets Hu and Hd, which needs to be properly 
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fine-tuned to have a small mass term mj^: 

H = -cos/3eH^ + smf3Hu. (2) 

At the high scale Ms the boundary conditions on the quartic Higgs coupling and on the Higgs- 
higgsino-gaugino couplings of the effective theory are determined by SUSY invariance and yield: 

X{Ms) = ^ [g\Ms) + g'\Ms)] cos^ 2/3, (3) 

guiMs) =g{Ms) sin /3 , gd{Ms) = g{Ms) cos /3 , (4) 

~g:,{Ms)=g' (Ms) sin P , g^iiMs) = g' (Ms) cos P , (5) 

where g and g' are the SU{2)l and U{1)y gauge couplings, tan/? is interpreted as the angle that 
rotates the two Higgs doublets into one heavy and one light (SM-like). 

Since for the MSSM with heavy scalars the number of basic input parameters is rather small, 
one can relax the assumption of a universal gaugino mass at Mqut and still have a predictive 
model which could lead to a different phenomenology with respect to the universal scenario. In 
reference [20], various scenarios for non-universal gaugino masses have been discussed. A particular 
interesting one is the gravity-mediated SUSY-breaking scenario in which the gaugino masses arise 
from a dimension-5 operator C oc {F^)ab/Mpi ■ A'^A'', where A are the gaugino fields and F$ is the 
auxiliary component of a left-handed chiral superfield $ which couples to the SUSY field strength. 
In the usual minimal Supergravity (mSUGRA) model with SU{5) grand unification, the field F$ 
is a singlet under the unifying gauge group, leading to universal gaugino masses. However, the 
field $ could sit in any representation of the symmetric product of the adjoint group [28j. For 
SU{5) symmetry, F$ could belong to an irreducible representation which results from the symmetric 
product of two adjoints 

(24 24)^ymmetric = 1 © 24 © 75 200 . (6) 

Once the neutral component of has acquired a vev, {F^)ab = Va^ab, the vevs Va determine 
the relative magnitude of the soft SUSY-breaking gaugino mass parameters Ma at Mqut [29l [30] . 
This is shown in the left-hand side of table [1] and, as can be seen, only in the singlet case 1 these 
parameters are universal. We also give the Ma values at the electroweak scale Mz, by using the 
one-loop renormalization group equations. In this vein, for a given representation besides the SM 
basic input parameters, the free parameters of the model are the common scalar mass Ms, tan/3, 
niMz) and Mi; M2 and M3 being fixed at Mqut- 

For the determination of the mass spectrum, the resummation to all orders of the leading 
logarithms of the large scalar mass and the one-loop radiative corrections has been taken into 
account, using the routine SHeavy of the code SuSpect [31] described in reference [20j. 

2.2 Dark matter within the MSSM with heavy scalars 

As deduced from the WMAP satellite measurement of the temperature anisotropics in the Cosmic 
Microwave Background, cold dark matter makes up approximately 23% of the energy of the Universe 
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Q = Mgvt 


Q = Mz 


1 


1:1:1 


1.0 : 2.0 : 6.2 


24 


1 : 3 : -2 


1.0 : 6.0 : -12.7 


75 


5 : -3 : -1 


1.0 : -1.2 : -1.3 


200 


10 : 2 : 1 


2.5 : 1.0 : 1.6 



Table 1: The ratios of gaugino mass parameters, Mi : M2 :M3, at the renormahzation scales Mgut 
and Mz, for the different patterns of soft SUSY breaking. The values Ms = 10^ GeV and tan /? = 5 
has been used. 



The DM cosmological density is precisely measured to be 

0£,A^/i2 = 0.1099 ±0.062, (7) 

which leads to 0.104 < Vt^Mh? < 0.116 at 68% CL. The accuracy is expected to be improved to 
the percent level by future measurements at Planck satellite |33| . 

As it is well known, the LSP neutralino can be an ideal WIMP-like cold DM candidate and in 
some areas of the SUSY parameter space the cosmological relic density associated to the lightest 
neutralino which is inversely proportional to the neutralino annihilation cross-section cJann = 
f^CXiXi ~^ SM particles), falls into the range required by WMAP. In the MSSM with heavy scalars, 
the neutralino mass matrix in the bino, wino, higgsino basis {B, W3, iJ^, H^) reads as follows: 
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where v ~ 174 GeV is the SM-like Higgs boson vev. Thus the lightest neutralino x? is a linear 
superposition of the bino, wino and higgsinos states: 
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NuB + NuWs + NisH^^ + NuH^; (9) 

the matrix rotates the neutralino states so that the mass matrix N Mn A^"^ is diagonal. 

In the MSSM with heavy scalars, there are essentially only four regions in which the WMAP 
constraint is fulfilled: 

- The 'mixed region' in which the LSP is a higgsino-gaugino mixture [34] , Mi ~ \fi\, which 
enhances (but not too much) its annihilation cross-sections into final states containing gauge 
and/or Higgs bosons and top quarks, XiXi ~^ W~^W~, ZZ, HZ, HH and ti. 

- The 'pure higgsino' and 'pure wino' regions, in which the LSP is almost degenerate in mass 
with the lightest chargino (xf) and the next-to-lightest neutralino (X2)- Such a scenario 
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leads to an enhanced annihilation of sparticles since the Xi ~ Xi and Xi — X2 coannihilation 
cross-sections |35] are much larger than that of the LSP. This solution generally requires LSP 
masses beyond 1 TeV. 

- The 'ff-pole' region in which the LSP is rather light, m^o ~ ^Mu, and the s-channel H 
exchange is nearly resonant allowing the neutralinos to annihilate efficiently |36l [37] . 

- The 'Z-pole' region in which the LSP is very light, m^o ~ ^Mz ~ 45 GeV, and the s-channel 
Z exchange is nearly resonant. We note that this region is not ruled out only in scenarios 
where the mass splitting between Mi and M2 at the electroweak scale is very large [20] , 

The computation of dark matter relic density has been performed using an adapted version of 
the micrOMEGAs public code [38l[39], where all the scalar particles (except for the SM-like Higgs 
boson) have been integrated out. Figure [H displays the area, in the [Mi, /i] plane, in which the 
WMAP constraint is satisfied (red-dark gray); a common scalar mass M5 = 10^ GeV, tan/3 = 5 
and the universal scenario has been chosen. Hereafter Mi and fj, have to be interpreted as parameters 
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100 1000 

Mj [GeV] 

Figure 1: The regions of the [Mi, fj] plane in which the WMAP constraint is fulfilled (red-dark 
gray area) for a common scalar mass value M5 = 10^ GeV, tan/? = 5 and the universal scenario. 
The green (light gray) area on the left and the bottom is the one excluded by direct searches of 
SUSY particles. 

evaluated at the electroweak scale. The green (light grey) area in the left and bottom parts of the 
figure denotes the region excluded by the collider data. This region is ruled out by the negative 
search of charginos at LEP2 from pair production of the lightest chargino: e~^e~ xfxi [40j. The 
peak for small Mi values at Mi ~ ^Mh, is due to the s-channel exchange of an almost real Higgs 
boson, Xi Xi ~^ H . For the mass value obtained here, Mh ~ 130 GeV, the Higgs boson mainly 
decays into hh final states. Right into the peak one is too close to the Higgs mass pole, and the 
LSP annihilation is too efficient leading to a too small VL^m h?- The peak reaches up to /i ~ 1300 
GeV, a value beyond which the LSP is almost bino-like and its coupling to the Higgs boson is too 
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small (the Higgs prefers to couple to a higgsino-gaugino mixture) to generate a sizable annihilation 
cross-section. 

For larger ^ and Mi values there is an almost straight band in which /j. ~ Mi ~ and 
the LSP is a bino-higgsino mixture with sizable couplings to W, Z and Higgs bosons, allowing for 
reasonably large rates for neutralino annihilation into 

X?X? ^ ZZ, HZ and HH final 

states. For instance, for Mi ~ 200 GeV and /i ~ 200 GeV, the annihilation cross-section is mostly 
due to the WW and HH final states and, to a lesser extent, the ZZ and ZH final states. For 
slightly larger /x and Mi values there is a jump due to the opening of the XiXi ~^ channel, which 
henceforth dominates the annihilation cross-section. Above and below the band, the LSP couplings 
to the various final states are either too strong or too weak to generate the relevant relic density. 
For fi values close to 1 TeV and even larger values of Mi there is a wider area in which the WMAP 
constraint is also fulfilled. In this region the LSP is almost a pure higgsino and a correct JIdm^^ 
can also be obtained thanks to the coannihilation of the LSP with the xf and X2 states. For lower 
fj, values and Mi still very large the LSP coannihilation with xf and X2 exceedingly strong and 
leads to a too small Qdm h'^- 

In figure[2]we present a similar analysis to the previous one, with the only proviso that tan /? = 30 
(left pane) and Ms = 10^° GeV (right pane). The fi gure in the left-hand side shows similar features 




100 1000 100 1000 

Ml [GeV] Ml [GeV] 

Figure 2: The same as in figure [1] but for tan/3 = 30 (left pane) or Ms = W^^ GeV (right pane). 

to figure [H a worth noting difference is that the peak which is due to the s-channel Higgs boson 
exchange reaches up to /x ~ 600 GeV, owing to the fact that for small tan P values the LSP becomes 
bino-like faster than in the high tan f3 case and its couplings to the Higgs boson are thus smaller. 
On the other hand, for M5 = 10^*^ GeV, the Higgs peak is shifted to a slightly higher Mi value. 
Ml ~ ^Mh ~ 80 GeV. In this case the annihilation channel XiXi ^H^WW*^Wff gives a 
significant contribution to the total cross-section [SO]- Note that variations over Ms and tan/3 are 
primarily reflected in the Higgs peak, whereas the mixed region is almost insensitive. 

Figure [3] is similar to figure [U with the difference that we take two situations presenting non- 
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universality in the gaugino masses at Mgut: scenario 24 (left pane) and scenario 75 (right pane). 
The case considered in the left-hand side shows similar features to the one with the universal setup 
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Figure 3: The same as in figured] but for scenarios 24 (left pane) and 75 (right pane). 

explored in the foregoing, except that a new peak at Mi ~ ~ 45 GeV appears, as a consequence 
of the s-channel Z boson exchange. In this scenario, such a resonance is no longer excluded because 
of the large gap between Mi and M2 (M2 ~ 6.3 • Mi at the electroweak scale). In this case, the 
invisible decay of the Z boson, Z XiXii and the chargino production measurements exclude a 
corner in the parameter space corresponding to low values for Mi and fi [20| . 

In the scenario 75 one has Mi : M2 : M3 ~ 1 : —1.2 : —1.5 at the weak scale, so that the LSP 
is in general close in mass to the lightest chargino and the next-to-lightest neutralino, and hence 
coannihilation of these states plays a very important role. There is a thin line up to Mi < 1 TeV 
where the LSP is dominantly bino-like and the WMAP Qdm h'^ range is obtained by an efficient 
annihilation and coannihilation of x?, X2 xf- Moreover, another region giving rise to the correct 
relic density corresponds to a large band with ~ 1 TeV and Mi > 1 TeV, where the LSP is almost 
a pure higgsino state. Once again, coannihilation leads the dark matter relic density. 
In scenario 200 the regions fulfilling the dark mater relic density compatible with WMAP constraints 
correspond to a LSP with mass > 1 TeV [20j. Owing to the fact that such a heavy candidate is 
largely outside capabilities of the next generation dark matter direct detection experiments, such 
scenario will not be considered in the analysis presented in sections H] and O 

3 Direct detection 
3.1 DifFerential event rate 

In spite of the experimental challenges, a number of efforts worldwide are actively pursuing to 
directly detect WIMPs with a variety of targets and approaches. Many direct dark matter detection 
experiments are now either operating or in preparation. All these experiments measure the number 
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N of elastic collisions between WIMPs and target nuclei in a detector, per unit detector mass and 
per unit of time, as a function of the nuclear recoil energy Er. The detection rate in a detector 
depends on the density po ~ 0.3 GeV cm~^ and velocity distribution f{vy-) of WIMPs near the 
Earth. Usually, the motion of the Earth is neglected and a Maxwellian halo for WIMP's velocity is 
assumed. In general, the differential event rate per unit detector mass and per unit of time can be 
written as: 



dN _ a^.N Po ^ ,2 r /(%) 



dv^ , (10) 



where the WIMP-nucleus scattering cross-section, ct^^n, is related to the WIMP-nucleon cross- 
section, <J^-p, by cr^^-Tv = <^x-p ' i^iT^r/Mr)'^ , with Mr = the WIMP-nucleon reduced mass, 
rrir = ^^]^^ the WIMP-nucleus reduced mass, the WIMP mass, niN the nucleus mass, and A 
the atomic weight. F is the nuclear form factor; in the following analysis the Woods-Saxon Form 
factor will be used. The integration over velocities is limited to those which c an give place to a 
recoil energy E^., thus there is a minimal velocity given by Vmin{Er) = 

In order to compare the theoretical signal with the background it is necessary to calculate the 
Let us call iV^*s" the signal, N'^'^^ the background and iV*°* = Ar^*9" + iv6fc9 the total signal measured 
by the detector. We will divide the energy range between 4 and 30 keV in n = 7 equidistant energy 
bins. For the discrimination between the signal and the background we calculate the variance y^: 



^[^Vr-Af^^j _ ^^^^ 



Here we are assuming a Gaussian error ai = \J jpj^ on the measurement, where M is the detector 
mass and T the exposure time. We use a similar analysis to that performed in references [24ll4T| . 

3.2 A Xenon-like experiment 

The Xenon experiment aims at the direct detection of dark matter via its elastic scattering off xenon 
nuclei. It allows the simultaneous measurement of direct scintillation in the liquid and of ionization, 
via proportional scintillation in the gas. In this way. Xenon discriminates signal from background 
for a nuclear recoil energy as small as 4.5 keV. Currently a 10 kg detector is being used, but the 
final mass will be 1 ton of liquid xenon. In figure [H we show the sensitivity curve for XenonlO 
(M = 10 kg) and XenonlT {M = 1 ton) for T = 3 years of data acquisition. 

In our study, following reference [42] we will always consider 7 energy bins between 4 and 30 
keV and 3 years of data acquisition for a 100 kg Xenon experiment. We could take into account 
non-zero background using simulations of the recoil spectra of neutrons in our analysis, and this 
would significantly degrade the sensitivity of the detector. However, this would involve a much 
more detailed study of the detector components (shielding, etc.), and we will not carry it out. In 
that sense, our results will be the most optimistic ones. Comprehensive studies about the influence 
of astrophysical and background assumptions can be found in references [241 ITOj . 
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Figure [5] shows the abihty of Xenon to determine the mass and scattering cross-section for a 
20, 100 and 500 GeV WIMP with cross-sections of 10~^ and 10~^ pb, in a microscopically model 
independent approach. We can clearly see how sensitive the experiment is to light WIMPs: the 




10 100 1000 

[GeV] 

Figure 5: Distribution of the maximum in the [o'x-p) it^x\ pl^ne, for 3 years of exposure in a 100 
kg Xenon experiment, for rriy^ = 20, 100, 500 GeV and cr-^-p = 10~®, 10~^ pb. The lines represent 
the 68% CL region, and the crosses denote the theoretical input parameters. 

precision can reach the percent level for < 50 GeV and high cr^-p. For WIMPs much heavier 
than the nucleus mass (~ 100 GeV for Xenon), the differential rate in equation (fTUl) becomes almost 
independent on and therefore the relative errors enlarge. However, for smaller values of the 
scattering cross-section < 10~^, the possibility of reconstructing or cr^-p, at least in a model 
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independent framework, quickly vanishes because of the dramatic increase of the errors. One way 
to deal with this deterrent is e.g. to place ourselves into the framework of a given model so as to 
limit the phase space (i.e. the {m^ , cFx-p) pairs) upon the inclusion of the available theoretical and 
experimental constraints. In this vein, hereafter we leave aside the model independent approach to 
the dark matter direct detection in order to focus to the MSSM with heavy scalars framework. 

4 Direct detection prospects within the MSSM with heavy scalars 

As far as the DM direct detection concerns, the cross-section for elastic scattering of a WIMP with 
a nuclei detector is perhaps one of the most important properties. This cross-section determines the 
detection rate in direct-detection experiments. The WIMP-nucleus elastic scattering cross-section 
depends fundamentally on the WIMP-quark interaction strength. 

In the MSSM with heavy scalars, the leading processes which give rise to the neutralino-nucleus spin- 
independent interaction are shown in figure [6l The diagrams correspond to the elastic scattering 
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Figure 6: Feynman diagrams contributing to the spin-independent elastic scattering of neutralinos 
with quarks and gluons, within the MSSM with heavy scalars. 

of a LSP from a quark and a gluon by the exchange of a Higgs boson in the t- and u-channels. 
The Higgs boson-gluon interaction is induced at the quantum level primarily by top-quark loops. 
Since both of the diagrams depend on the interaction between the lightest neutralino and the Higgs 
boson, the XXH vertex will drive the phenomenology of the scattering cross-section. This coupling 
is given by: 

C^o^oH Of ^^13 {9d Nu - 9d Nu) - Nu {~9u Nu - ~gi iVii) • (12) 

As equation (fT2]) depends on the medley of the whole Nij matrix elements, it follows that the Hxx 
coupling can be enhanced for a lightest neutralino which is a 'temperate' gaugino-higgsino mixture. 
In addition, for a pure higgsino-like or a pure gaugino-like neutralino, the coupling with a Higgs 
boson vanishes, and hence it does the scattering cross-section. 

In the general MSSM case, there also exists a potentially important input to the cross-section 
coming from a plethora of diagrams involving squarks propagators (essentially the first generation 
squarks u and d), both contributing to the interaction of WIMPs with quarks and with gluons [l]. 
These processes could be dominant in the case when the squarks are light. Nevertheless, in this case 
we are taking into account heavy scalars, hence diagrams containing squarks will be suppressed by 
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powers of 

In the same vein, let us recall that the Z-boson exchange in the interaction XiQ ~^ XiQ does not 
contribute to the spin-independent scattering cross-section but only to the spin-dependent one. As 
already pointed out [13], spin-dependent direct neutralino searches in the general MSSM appear to 
be largely disfavored with respect to spin-independent ones. 

In the same way as for the dark matter relic density, the computation of scattering cross-sections 
has been performed using an adapted version of the micrOMEGAs code. Thus, it is possible to 
use equation (fTOl) in order to extract the expected differential event rate of nuclear recoil, in the 
framework of the Xenon experiment. Figure [7] presents the contour lines for the scattering cross- 
section cr^-p (dashed lines) and the mass of the lightest neutralino (dotted lines) in the [Mi, ^] 
plane, for Ms = 10^ GeV, tan (3 = 5 and the universal scenario. It is also shown the region allowed 




100 1000 
Mj [GeV] 



Figure 7: Contour lines for the scattering cross-section cr^^p (black dashed lines) and the mass of 
the LSP (blue dotted lines) in the [Mi, /x] plane, for Ms = 10^ GeV, tan/? = 5 and the universal 
scenario. The regions fulfilling the WMAP constraint (red-dark gray) and excluded by collider 
searches (green-light gray) are also shown. 

by the WMAP measurements and the parameter space excluded by collider data. Worth noticing 
is the fact that while Mi ~ fi, the cross-section reaches high values, up to ~ 4 • 10~^ pb. The 
XXH coupling stays high even for elevated Mi and fj, values, but the maximum slowly narrows. 
However, as we move away from the higgsino-gaugino mixing region, <Tx-p decreases severely. Let 
us emphasize that in the Higgs-pole region, the bino-like nature of the lightest neutralino implies a 
highly suppressed scalar cross-section via Higgs boson exchange because of the low xxH coupling. 
In that sense, the Higgs boson resonance does not enhance the scattering cross-section (T^-p- 
Changing Ms to higher values causes a rapid decrease in the scattering cross-section. This may be 
largely explained by the variation of the Higgs boson mass, which increases logarithmically with 
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Ms, and by the strong dependence of the cross-section c-^-p on the mass Mh'- 

^x-P ^ ■ (13) 

On the other hand, the increase of tan /? reduces the scattering cross-section because of two effects: 
i) the slight growth of the Higgs mass and ii) the LSP becomes more quickly a pure bino- or higgsino- 
like state. Smaller values for tan /? (< 2 — 3) produce a SM-like Higgs boson which becomes too 
light, actually excluded by negative searches at LEP2. In this vein, the analysis in the next section 
will be focused to the kind of best scenario for dark matter direct detection, i.e. Ms = 10^ GeV 
and tan/3 = 5. 

The maximum scattering cross-section reached in figure [7] is not yet ruled out by dark matter 
direct detection experiments; however, some collaborations aim to exclude cross-sections up to lO"^'^ 
pb in a not too distant future. In figure [H] we show exclusion lines for Xenon after 3 years of data 
acquisition, on the [Mi, /j] parameter space, in the framework of the universal scenario. These 
curves reflect the Xenon sensitivity and represent its ability to test and exclude different regions of 
the MSSM with heavy scalars at 68%, 95% and 99% CL. The regions which could be excluded follow 
approximately the contour lines for the scattering cross-section (see figure [7]). Moreover, they tend 
to enlarge for LSP masses of the order of m-^ ~ 100 GeV, as expected from figure [H As can be seen 
in figure m for the universal scenario the absence of signal at Xenon could exclude a sizeable fraction 
of the parameter space of the MSSM with heavy scalars, in particular a major part of the region 
fulfilling the WMAP constraint. In fact, only the top of the Higgs peak and the area corresponding 
to ~ 1 TeV and Mi > 1 TeV could not be probed. This is due to the low scattering cross-section 
and the high LSP mass respectively. 

Figure M is similar to the previous one, with the difference that we explore the non-universal 
scenarios 24 and 75. In the same way as for the universal case, for scenario 24 the exclusion lines 
follow approximately the /i ~ Mi line, getting large for WIMP masses of the order of 100 GeV. For 
"^x ~ ^0 GeV, the sensitivity of the experiment deteriorates rapidly (see figure HJ. Let us note that 
the dark matter scattering cross-section is not enhanced by the resonant exchange of a Z boson, 
in the same way as for the Higgs boson peak. Finally, the direct detection prospects for scenario 
75 are not encouraging. In fact, in this case the regions leading to the correct dark matter relic 
density measured by WMAP are not generated by the usual annihilation process, but rather by 
coannihilation processes of the LSP with the X2 ^^'^ the xf- Hence, most of the parameter space 
region that a Xenon-like experiment could examine corresponds to an annihilation cross-section 
which turns out to be too high indeed to be measurable. Nevertheless, a small single region exists 
with fi ~ Ml ~ 1 TeV which gives rise to the correct relic density and which could be potentially 
excluded by a Xenon-like experiment; but such region corresponds to a very heavy LSP, at the 
actual border of detectability. 



12 




100 1000 100 1000 

Ml [GeV] Ml [GeV] 



Figure 8: Exclusion lines for 3 years of exposure in a 100 kg Xenon experiment, on the [Mi, ^] 
plane, for different combinations of Ms and tan/3, in the framework of the universal scenario. The 
dashed lines represent the 68%, 95% and 99% CL exclusion lines. The regions fulfilling the WMAP 
constraint (red-dark gray) and excluded by collider searches (green-light gray) are also shown. 

5 Reconstruction prospects 

Once a signal has been detected in a given direct dark matter detection experiment, a new question 
appears concerning the possibilities of reconstructing both physical observables and fundamental 
parameters of the model. Gradually, a new generation of dark matter experiments start to make 
measurements and not simply to set limits [Sj. In this vein, the aim of this section is to assess, 
first of all, the possibilities of identifying neutralino LSP properties, as its mass and its scattering 
cross-section, in the framework of the MSSM with heavy scalars from a Xenon-like experiment. 
Secondly, we will also examine the prospects of reconstructing the fundamental parameter space, 
i.e. the low-energy Lagrangian parameters. Even if in the general MSSM it is a formidable task, the 
reduced number of free parameters of the present model makes the work more feasible. However, for 
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reconstructing all parameters using direct detection experiments there is a first difficulty because 
they only give rise to a single observable: the differential event rate of nuclear recoil. This can 
be partially eased by taking into account the existing constraints coming from collider physics but 
especially the one given by WMAP. A second complication comes from the fact that the sector of 
the parameter space that can be tested by Xenon is almost insensitive to some of the parameters. 
Actually, Ms and tan/3 will essentially determine the mass of the SM-Higgs boson, but their impact 
in the region where the LSP is a higgsino-gaugino mixture is somehow limited (e.g. see figures 
[UI2] and [8]) . Therefore, we will concentrate on the reconstruction over the [Mi, /i] plane in the 
framework of a given pattern of soft SUSY-breaking gaugino masses. Other measurements (e.g. at 
LHC or ILC) could bring complementary information about the remaining low-energy Lagrangian 
parameters. 

In order to study the abilities of Xenon to reconstruct the parameter space of the MSSM with 
heavy scalars, we have chosen four characteristic benchmark points defined in table[2l corresponding 
to different values for the (Mi, /x) pair and different scenarios for gaugino masses. We have stuck to 





Scenario 


Ml [GeV] 


/i [GeV] 


mx [GeV] 




P [Pb] 


A 


1 


138 


143 


93.6 


3.2 


•10-9 


B 


1 


207 


262 


185.1 


1.7 


•10-9 


C 


1 


71 


243 


63.2 


3.0 


10-10 


D 


24 


45 


165 


39.0 


2.9 


10-10 



Table 2: Benchmark points used throughout the analysis, for M5 = 10^ GeV, tan/3 = 5 and 
different scenarios of gaugino masses. The corresponding LSP mass and scattering cross-section are 
also presented. 

low values for the parameters Ms and tan/3 (M5 = lO'* GeV and tan/? = 5) in order to maximize 
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the direct detection potential of Xenon. Let us recall that in this limit the scattering cross-section 
is maximized. All these points satisfy the whole collider constraints and also generate a dark matter 
relic density in agreement with WMAP measurements. Moreover, they belong to the region that 
could be probed by a 100 kg Xenon experiment after 3 years of exposure. For each point, table [2] 
also presents the corresponding LSP mass and scattering cross-section. 

Points A and B correspond to the mixed region in which the lightest neutralino is a higgsino-bino 
mixture and then, they give rise to a high scattering cross-section cr^^p > 10"^ pb. However, from 
the point of view of direct detection, point A presents an important advantage over point B. In fact, 
point A generates a neutralino with mass rriy. ~ 100 GeV, corresponding to the maximal sensitivity, 
while point B produces a heavier LSP. In the latter case, the dark matter relic density is mainly 
produced by the reaction Xi Xi ^ ^ t- On the other hand, point C corresponds to the Higgs pole 
region in which the lightest neutralino is rather light, ~ ^Mh ~ 65 GeV. Despite its low mass, 
this point gives rise to a very weak scattering cross-section cr^-p ~ 10"^'' pb, in the threshold of 
detectability. Finally, point D is located over the Z-pole peak and corresponds to an LSP with a 
very light mass (m^ ~ ~ 45 GeV) and scattering cross-section. Because of the large mass 

splitting between the bino and the wino mass terms in scenario 24, M2 ~ 6.3 -Mi at the electroweak 
scale, such an LSP mass is not excluded. 

Figure [10] shows the ability of Xenon to reconstruct the Mi and /x parameters, for the four 
benchmarks described in table El Even if the reconstructed regions are in general relatively large, 
in particular for points B and C, the conjunction with collider and cosmological constraints allows 
to drastically shrink the latter. 

For point A, for instance, the overlap of the different bounds strongly restricts the region that can- 
not be discriminated by Xenon in the [Mi, jj] plane to Mi ~ 141 GeV and 121 GeV < /i ^ 165 GeV. 
Henceforth, the regions will correspond to a 68% CL. This range corresponds to a lightest neutralino 
with a mass within the interval 84 GeV < < 104 GeV and a cross-section cr^^p between 2.9-10"^ 
pb and 3.4 • 10~^ pb. Note that, for this case, both the LSP mass and the scattering cross-section 
can be very well reconstructed with a precision of the order of 20% and 15% respectively. 
With regard to points B and C, they present highest values for the LSP mass or the cross-section 
implying a spoilage of the reconstruction capacities. Actually, for those points, the discrimination 
regions grow and branch in two areas almost symmetrical about the /i = Mi axis. The recon- 
structed area compatible with constraints for point B consists of two parts. The largest one almost 
follows the fj, = Ml line with Mi varying within the interval delimited by 197 GeV and 377 GeV 
corresponding to a LSP with a well defined scattering cross-section of the order of ~ 2 • 10"^ pb 
but with a mass within the interval of 173 — 350 GeV. Moreover, the signal is also compatible with 
a very narrow region located over the Higgs-peak and corresponding to Mi ~ 75 GeV and 172 GeV 
< /X < 192 GeV. This region gives rise to a ~ 61 GeV LSP with a scattering cross-section between 
1.1 • 10~^ pb and 7.5 • 10~^° pb. A 100 kg Xenon experiment after 3 years of exposure would not 
be able to disentangle these two unrelated areas. 

The region engendered by point C is located around Mi ~ 71 GeV and 218 GeV < /x ^ 285 GeV. 
In such an interval the mass of the lightest neutralino is close to 61 GeV but the cross-section 
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100 1000 100 1000 

Ml [GeV] Mj [GeV] 

Figure 10: Ability of a 100 kg Xenon experiment after 3 years of exposure to reconstruct the Mi 
and fi parameters, for the three benchmarks of the tabled at 68% and 95% CL. The crosses denote 
the theoretical input parameters. The regions fulfilling the WMAP constraint (red-dark gray) and 
excluded by collider searches (green-light gray) are also shown. 



can vary in a large range, from 1.5 • 10~ pb to 4.6 • 10~ pb. There is, nevertheless, an extra 
region compatible both with the dark matter signal and the whole set of constraints, standing near 
Ml ~ |i ~ 0.7 — 1 TeV. However, even if this area gives rise to a large scattering cross-section of 
~ 2 • 10~^ pb, because of the large mass of the LSP, we are in a regime where the detector is at the 
limit of its capabilities. 

Finally, the reconstructed region fulfilling the constraints for point D consists of four parts corre- 
sponding to the left and right bands of the Z- and ff-peaks. They are roughly limited by 160 GeV 

< /i < 255 GeV and Mi ~ 45 GeV, 55-60 GeV and 68 GeV. In the four areas, the LSP mass varies 
in the range between 37 GeV and 62 GeV; the scattering cross-section varies between 1.2 • 10"^'' pb 

< c^x-p ^ 3.4 • lO"^'^ pb. For this case, the LSP mass can be reconstructed with a precision of the 
order of 65%, whereas, for the cross-section the relative error reaches almost 75%. It is clear that 
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a Xenon-like experiment cannot examine with a high-precision level such a scenario with a so low 
scattering cross section. However, it can provide very valuables hints on the nature of the WIMP 
dark matter. 



6 Conclusions 

We have explored the dark matter detection prospects in the Minimal Supersymmetric Standard 
Model in the scenario where the scalar partners of the fermions and the Higgs particles (except for 
the Standard-Model-like one) are assumed to be very heavy and are removed from the low-energy 
spectrum. In the MSSM with heavy scalars, the WIMP candidate for the dark matter relic density 
is the neutralino LSP. We have analysed the neutralino LSP (x?) in scenarios where the gaugino 
mass parameters are universal but also the case where they are non-universal at the GUT scale. 
This analysis has been carried out in the framework of a Xenon-like 100 kg experiment. In general, 
an important fraction of the parameter space giving rise to the dark matter relic density measured 
by WMAP can be probed and excluded in the case of not detecting any WIMP. In the opposite 
case, once a WIMP signal has been found, we have shown that for a light Xi which is a higgsino- 
gaugino mixture it is possible to reconstruct efficiently the mass and the scattering cross-section of 
the neutralino LSP. Moreover, we have shown that it is also feasible to put strong constraints over 
some of the parameters of the Lagrangian, e.g. the higgsino and the gaugino mass parameters. 
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